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Abstract 

The behavior of the geometric phase gained by a single spin- 1/2 nucleus immersed into a thermal 
or a squeezed environment is investigated. Both the time dependence of the phase and its value at 
infinity are examined against several physical parameters. It is observed that for some intermediate 
ranges of the temperature and the coupling strength, the presence of squeezing enhances the 
geometric phase. 
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I. INTRODUCTION 



The geometry of the Hilbert space of a quantum system is registered to the memory of the 
system as a geometric phase factor [1, 2, 3, 4, 5, 6, 7]. It is inherently fault tolerant because of 
its geometric nature [8]. For this reason, it is used in the implementation of controlled gates 
for quantum computers[9, 10]. These gates are also of crucial importance for a universal 
set of quantum logic gates[9, 10, 11]. For this reason, a full understanding of the nature of 
quantum geometric phase in different environments is vital for quantum computations. 

Various aspects of the effects of the environment on the GP of open quantum systems 
have been studied. Rezakhani and Zanardi analyzed the temperature effects on mixed-state 
GP for a single and two coupled spin-1/2 particles[12]. Wang et al. analyzed the effects 
of a squeezed vacuum reservoir on GP of a two-level atom in an electromagnetic field by 
a formulation entirely in terms of geometric structures [13]. Carollo et al. showed that 
GP can be induced by cyclic evolution in an adiabatically manipulated squeezed vacuum 
reservoir [14]. Banerjee and Srikanth studied the effects of a squeezed-thermal environment 
on the GP of a two-level system[15] for dissipative and non-dissipative cases, and analyzed 
the initial-state and temperature dependence of GP for the system. The purpose of this 
article is to examine the effect of various physical parameters on the time-dependence of 
GP. 

The system analyzed in this article is a two-level nucleus inside a magnetic field. The 
time evolution is non- unitary due to the interaction with the environment. The environment 
is first taken as a thermal bath. Then, this bath is considered to be driven by an electro- 
magnetic field in a squeezed state in order to see whether the GP can be enhanced. The 
dependence of the GP on temperature, external magnetic field, coupling strength, squeez- 
ing and initial state of the system are then analyzed. The GP is computed by using the 
kinematic definition given by Tong et al. [16]. In their approach, when the system's density 
matrix in the Schrodinger picture is pit), the GP gained in the time interval [0, r] is given 
by the expression 

$(r) = arg ^ v / A fc (0)A fe (r)(^(0)|^(r))e(-^(^WI^W)^) J , (1) 

where \k(t) and | are the eigenvalues and the eigenvectors of pit), respectively. 

The content of the article is as follows. In section II, the system and its interaction with 
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the environment is described. This section also includes the derivation of the non-unitary 
dynamics of the reduced density matrix by using Markov approximation. After that, the 
density matrix of the system is computed analytically. In section III, the dependence of the 
GP on various physical parameters is analyzed. Finally, section IV contains brief conclusions. 



II. THE NUCLEUS INSIDE A BATH 

Our specific system is a single spin-1/2 nucleus in an external static magnetic field B 
which is taken to be in the z direction. The Hamiltonian of the nucleus is 

H N = -huj N i z , (2) 

where h is the Planck constant, ujn = 7at|-B|, 7^ is the gyromagnetic ratio of the nucleus, 
I z = a z /2 and a z is the z— component of the Pauli spin operator. 

The environment (reservoir) is assumed to be a bath of harmonic oscillators (such as 
electromagnetic radiation) where the frequency spectrum forms a continuum. The annihi- 
lation (creation) operator for the mode at frequency ui is denoted by b{uj) {W{uj)) and the 
Hamiltonian of the environment is 

poo 

H R = h ub\u)b(u)6u . (3) 
Jo 

The interaction of nucleus with the oscillators is described by the interaction Hamiltonian, 

POO 

H NR = h g N (uj)I + b ] (uj)duj + h.c. , (4) 
Jo 

where Qn{^) is the coupling coefficient between the nucleus and the cu-mode of the reservoir, 
= I x ±il y are the spin ladder operators, and h.c. indicates the hermitian conjugate term. 
The total Hamiltonian H = Hn + Hr + H^r-, is transformed into the interaction picture 

as 

poo 

V(t) = h / gN^e^^i+P^du +h.c. . (5) 
Jo 

Let f>Nn{t) denote the state of the combined system of the nucleus and the environment at 
time t, and let pi/n{t), p R {t) and p^{t) denote the interaction picture states of the combined 
system, the environment and the nucleus, respectively. We study this system in the regime 
in which the Markov approximation can be applied. In this approximation, it is assumed 
that the environment remains in its initial state during the evolution, p R (t) = Pr{0), and 
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the state of the combined system is taken to be in product form, plf R (t) = p^it) (£> pn(0). 
As a result, the equation of motion of the state of the nucleus can be obtained as 

p f N (t) = -^tr R [y(t)X(0) ®p R (0)} - ^rJ q [V(t), [V(t') rpW) ® Pfl(O)] ] < (6) 

where the dot denotes the time derivative and tr# represents trace over the degrees of freedom 
of the environment [17, 18]. 

For the thermal environment, the state is given by Pr{0) = pth = 
exp(—Hft/kBT)/trexp(—Hn/kBT)). If the bath is driven by a squeezed field, the state 
of the squeezed-thermal environment is given by 

Pr(0) = Sfa&, (7) 

where 

S = S[£(cu)} = exp j duo[i(uo)b\uj)b\2n - uj) - C(uj)b(2Q - ou)b(cu)]^j , (8) 

where 2Vt is the squeezing- carrier frequency, £,{uj) = r(oj)e i ^ u '\ r(uj) and 4>(uj) are real 
numbers characterizing the squeezing which satisfy = £(2f2 — uj). 

For both types of environments, we have (b(uj)) = which makes the first term on the 
right-hand side of Eq. (6) vanish. The second term can be expressed in the form 

~p N {t) = - l -[AH N ,p N (t)] + c_ (2i-p N (t)i + - p N (t)i + r - f + rp N (t)^ 
+c + (2f + p N (t)f- - p N (t)i~f + - i~i + p N {t)) 
-d i + p N (t)i + - d* tp N (t)r, (9) 

where AH^ — —hAuj^a z /2 represents a re-normalization of the frequency uj^ of the nucleus, 
and C± and D are numbers that capture the collective effect of all of the oscillators in the 
reservoir satisfying 

C + = C^ + it\g N {uj N )\ 2 . (10) 

Each of these constants can be expressed as an integral over the frequency uj, which can 
be decomposed into a principal-value integral, which depends on the precise frequency de- 
pendence of gN^), and a Dirac-delta integral term, which depends only on g N (uj N ). The 
coupling function ^(cu) is freely adjustable; by changing that function suitably, the values 
of all principal-value integrals and therefore all final constants C±, D and Aujn can be ad- 
justed to any desirable value. Here, for simplicity, all of these principal-value integrals are 
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taken to be zero and the resultant constants are used. With this choice, the re-normalized 
Hamiltonian AHn becomes zero. 

For the purely thermal environment, the constants in Eq. (9) can be computed as D = 
and C- = Tx\g^{uN)\ 2 n{uN), where 

< u) = exp(^/LT) - 1 

denotes the average occupation number of the mode at frequency u. In this case, the matrix 
elements of the interaction picture density matrix can be integrated to 

Pmiit) = P^n(oo) + (p^ 11 (0)-p^ 11 (oo))e- 2 ^ +c -) i , (12) 
ftw® = PmMe- {C++C - ]t , (13) 

where 

C 

Pmxi 00 ) = Patii(oo) = c * c ■ ( 14 ) 
For the squeezed-thermal environment given in Eq. (7), the constants can be found as 



C- = k\qn{un)\ |^(^Ar) + [n(oj N ) + n(2tt — u N ) + 1] sinh r(uj N )j , (15) 
D = 7rg N (u N )g N (2Q - u N ) [n(uj N ) + n(2fi - oo N ) + 1] sinh(2r(c; J v))e 2i(Q - a;iv) '-^ iv) (16) 

For this case too, the time-dependent density matrix can be found analytically. The diagonal 
entry p^r U (t) is still given by Eq. (12) and the off-diagonal entry can be expressed in general 
as 

pi m2 (t) = (A ie st + A 2 e~ st ) e -(c + +c-)t+i(a- Ulf )t ? (17) 
where s is the purely real or purely imaginary constant 



s = y/\D\*-(n-u N y, (is) 

and Ai and A 2 are constants that should be determined from the initial state Ptv(O). The 
squeezing changes the time-dependence of the density matrix as follows. First, it changes 
the long-time limit of the density matrix pjv(oo). It also makes the diagonal relaxation 
time, (2(C + + CJ))~ l , shorter. Apart from those, it changes the time-dependence of the 
off-diagonal entry; in particular it produces new oscillatory behavior for sufficiently large 
squeezing. 
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The geometric phase $(r) given in Eq. (1) can be computed analytically in the thermal 
case. The argument of the exponential in this expression is given by 

- j\ip±{t)\<P±{t))M = (t + £-1^-(F(t) - F(0))) (19) 



where 



F(t) = ln(A T -6 T + 25^) - In ( 2(5 - U + '^f^ ) , (20) 



A T = Ji-detpAr(r) , (21) 
$t = \ {pNii{r) - P N22{r)) , (22) 

and |v?±(t)) are eigenvectors with corresponding eigenvalues A±(r) = 1/2 ± A T . 

It should be noted that when r goes to infinity, the density- matrix Pn(j) goes to a 
diagonal time-independent state. Hence F{r) and consequently the phase $(r) tend to 
constant limits. For the squeezed-thermal environment, the argument of the exponential 
needs to be computed numerically but it can be shown that the same behavior holds at 
infinity, i.e., the GP settles down to a well-defined finite limit. 



III. THE DEPENDENCE OF THE GEOMETRIC PHASE ON PHYSICAL PA- 
RAMETERS 

For investigating the behavior of the GP under different physical conditions, the principal- 
value integrals are taken to be zero as explained in the previous section, and both of the 
coupling coefficients g N {u N ) and g N (2Q—u N ) are taken to be equal to a real positive value g. 
The phase parameter 4>(lon) of the squeezed state is taken to be 0. The relevant parameter 
r(uN) that gives the amount of squeezing is simply denoted by r below. 

In the problem, there are three parameters having the dimension of energy: the excitation 
energy hu^ (which is proportional to the external magnetic field), the thermal energy k^T 
and an energy related to the coupling strength hg 2 . The behavior of the GP is invariant if 
all three of these parameters are scaled by the same amount. For this reason, only the ratios 
of these parameters need to be specified as done below. 
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A. Time Dependence of the Geometric Phase 

The dependence of the geometric phase $(t) on time r has periodic structures with 
period equal to 2tt/uj^. In each of the oscillations over one period, 5>(t) also changes by an 
amount depending crucially on the values of all physical parameters. If there is squeezing, 
an additional periodic structure coming from the time-dependence of the off-diagonal entry 
given in Eq. (17) may appear. 

The effects of changing the temperature, the coupling strength and the magnetic field on 
the time-dependence of the GP are shown in the Figures 1, 2 and 3, respectively. As it can 
be seen, increasing the temperature or the coupling strength, or decreasing the magnetic 
field have regular effects when there is no squeezing (the figures on the left of Fig.s 1, 2 and 
3). With these changes, the oscillations of the GP are destroyed and the limiting value of 
GP is reached at earlier times. 

But when there is squeezing, that regular dependence is lost for the case of increasing 
the temperature or the coupling strength (Fig.s 1 and 2). For sufficiently large values of the 
temperature or the coupling constant, the limiting value of the GP is higher in comparison to 
those with lower temperatures or coupling strengths. Moreover the GP reaches the limiting 
values at a shorter time for higher temperatures and at a longer time for higher coupling 
strengths. For low temperatures and coupling strengths, presence of squeezing in general 
decreases the GP. 
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FIG. 1: The GP vs time for different temperatures for a thermal environment (left) and a squeezed- 
thermal environment with r = 1 (right). The coupling strength is g 2 /ujn = 0.01 and the initial 
state is a pure spin up state along x. The numbers on the graph indicate the value of fc^T '/hui^. 
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FIG. 2: The GP vs time for different coupling strengths for a thermal environment (left) and a 
squeezed-thermal environment with r = 2 (right). The magnetic field is hw^jhsT = 1 and the 
initial state is a pure spin up state along x. The values of the ratio g 2 /ujn are shown on the graphs. 
As can be seen, squeezing enhances the GP for a particular case with large coupling. Note also 
small oscillations which are produced by the effect of squeezed field on the spin state. 

For high magnetic fields, the squeezing has almost no effect (Fig. 3). At lower fields, the 
squeezing tends to decrease the GP, destroys the oscillations and make the limiting values 
reached at earlier times. And for small enough field strengths it is possible to get higher 
values for the GP with increasing the squeezing. 




FIG. 3: The GP vs time for different magnetic fields for a thermal environment (left) and a 
squeezed-thermal environment with r = 1 (right). The coupling coefficient is taken as kg 2 /ksT = 
0.01 and the initial state is a pure spin up state along x. The values of the ratio Hujn /^bT are 
shown on the graph. We have taken = 3lun for all data points. 
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Figure 4 shows the dependence of the GP on the initial state. The initial state is taken 
as a pure state along a direction on the xz-plane having the spherical angle 9. When 9 is 
increased above tt/2, the GP changes significantly. Especially, the GP tends to decrease at 
the initial moments for 9 > it/ 2. If the initial state is mixed, then the GP is suppressed in 
magnitude but its overall behavior does not change. 
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FIG. 4: The GP vs time for different initial states for a thermal environment (left) and a squeezed- 
thermal environment with r = 1 (right). The magnetic field is hio^lksT = 1 and the coupling 
strength is g 2 /ujn = 0.01. The initial state is chosen as a pure state and the value of the spherical 
angle 9 of the initial orientation is shown on the graph. 

B. The Dependence of the Phase at Infinity on Physical Parameters 

The behavior of the GP at infinity is also a quantity of interest. Its magnitude gives an 
indication of the overall magnitudes of the GP that can be obtained at any time. This value 
can also be used for showing the effect of the squeezing on the GP. In figure 5, the long time 
limit of the GP, i.e., $(oo), is shown as a function of temperature. When there is no squeez- 
ing, increasing the temperature always suppresses the GP. With the presence of squeezing, 
three different temperature regimes appear. For very low and very high temperatures, the 
GP is still suppressed with increasing T. However, there is now an intermediate range of 
temperatures over which the opposite trend appears. Increasing the squeezing parameter r 
moves this intermediate range to lower temperatures. Fig. 6 shows the dependence of the 
limiting value of the GP on the coupling strength. The behavior is similar to the case of 
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changing the temperature parameter. 
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FIG. 5: The GP at infinity as a function of temperature for different values of the squeezing 
parameter r. The coupling strength satisfies g 2 /u>N = 0.01. 
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FIG. 6: The GP at infinity as a function of the coupling strength for different values of the squeezing 
parameter r for the case fiWN jk^T = 1. 
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Figure 7 shows the dependence of the limiting value of the GP on the magnetic field. 
It can be seen that there are two regions for this case. For high enough magnetic fields 
this dependence is almost linear. Squeezing essentially decreases the slope. But at lower 
magnetic fields, GP has a different behavior under squeezing. In this regime, squeezing 
enhances the GP which is now a nonlinear function of the magnetic field. The limiting GP 
also reaches to a maximum value at an optimum value of the magnetic field. Increasing the 
squeezing expands that region of nonlinearity to higher magnetic fields. 
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FIG. 7: The GP at infinity as a function of the magnetic field for different values of the squeezing 
parameter r. The coupling strength satisfies hg 2 /ksT = 0.1. 

The dependence of $(oo) on the initial state is shown in Fig. Ill B. As can be seen, for 
low squeezing $(oo) reaches its largest values for 9 near tt, i.e., a spin orientation which is 
almost spin down. Small squeezing suppresses the GP but does not change the behavior. 
For large squeezing however, $(oo) can be enhanced significantly and its largest values are 
reached when 9 is around tt/2. 



11 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

0/tt 

FIG. 8: The GP at infinity as a function of the initial state angle 6 for different values of the 
squeezing parameter r. The parameters are chosen as g 2 /ujn = 0.025 and fkUN /ksT = 1. 

IV. CONCLUSION 

The availability of easy manipulation with the current technology NMR techniques makes 
nuclear systems good candidates as carrier systems of the entities that are necessary for 
quantum information processing[9, 19, 20, 21, 22, 23]. With this motivation, the effects of 
the coupling between a nucleus and a dissipative environment on the GP that the state of the 
nucleus gains during its time evolution are studied. Although it is affected by the physical 
parameters, it is always possible to get a finite GP. Its dependence on these parameters are 
almost as expected; it decreases with increasing temperature or coupling strength, because 
both increase the decoherence caused by the environment, and it increases with increasing 
magnetic field, since increasing B has the same effect as decreasing both the temperature 
and the coupling constant. 

If the GP needs to be increased further, the environment could be driven by an electro- 
magnetic field in a squeezed state. However, in this case, one should be very careful because 
the squeezing may have a destructive or an enhancing effect on the GP depending on the pa- 
rameters. In order to enhance the GP, the temperature should be held at some appropriate 
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intermediate value; it should not be too small or too large. The coupling strength should be 
taken at a relatively high value so that the squeezed field can interact with the spin more; 
but it should not be too large lest the spin relaxes to its equilibrium state quickly. And 
larger magnetic fields do not always enhance the GP; it is true only up to an upper bound 
for the field strength. A geometric phase that can be described by the model used in this 
article can promote the robustness in quantum computations. 
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